Energy Flow Puzzle of Soliton Ratchets 
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We study the mechanism of directed energy transport for soliton ratchets. The energy flow appears 
due to the progressive motion of a soliton (kink) which is an energy carrier. However, the energy 
current formed by internal system deformations (the total field momentum) is zero. We solve the 
underlying puzzle by showing that the energy flow is realized via an mhomogeneous energy exchange 
between the system and the external ac driving. Internal kink modes are unambiguously shown to 
be crucial for that transport process to take place. We also discuss effects of spatial discretization 
and combination of ac and dc external drivings. 
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A transport mechanism of potential relevance in var- 
ious areas of physics, chemistry and biology is based on 
the ratchet effect Q, i.e. the generation of directed cur- 
rents by zero-mean external perturbations. It is based 
on the breaking of relevant space-time symmetries of the 
underlying system evolution equations 0- A paradig- 
matic model corresponds to a classical particle moving 
in a spatially periodic potential under the influence of 
zero-mean fluctuations [l| , where the energy current is di- 
rectly connected with the mean particle momentum and 
the corresponding kinetic energy of the particle. For spa- 
tially extended systems, e.g. an annular Josephson junc- 
tion, which is described by a partial differential equation 
(PDE), the ratchet phenomenon manifests as a unidi- 
rectional motion of a collective kink excitation (soliton) 
0, 0, 0, 0, Hi ■ Here the unambiguous ab initio defini- 
tion of a current may become a much more complicated 
task, because the kink excites other modes in the system 
during its motion, which may contribute to an energy 
current as well. 

A well-known model in the field of soliton ratchets 
is the driven-damped sine-Gordon equation [^, which 
is also used for modelling the abovementioned annular 
Josephson junction j^: 



-a(ft — sinip + E{t) , 



(1) 



where E[t) is a zero-mean time-periodic driving force, 
E{t + T) = E{t), E{t)dt = 0. We impose the kink- 
bearing periodic boundary condition: 

ip{x + L,t) ^ ip{x, t) + Q , ipt{x + L,t) ^ ipt{x, t) , (2) 

where Q — 27rm is the topological charge with integer 
TO = 1, 2, and L is the system size. 

Let us consider the easiest case to = 1, i.e. the presence 
of one kink in the system Q — 2tt. The kink velocity V 
is defined e.g. as 0,13: 
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In the case of a soliton ratchet the mean value of V{t) will 
be nonzero. Since the kink carries some energy, one ex- 
pects a mean nonzero energy current as well. Recently it 
has been proposed to observe this directed energy trans- 
port using the definition of the internal energy current 
and its density Q : 



J\t) = / j^dx , j\x,t) = -ipxft ■ 



(4) 



is also known as the total momentum of the system 

miaiiii. 

Choosing either Eq.Q or Eq.Q), the symmetry anal- 
ysis provides identical necessary conditions for the ap- 
pearance of a ratchet effect 5]. If the ac driving E{t) 
possesses a shift symmetry. 



E{t) = -E{t + T/2) , 
then the combined symmetry transformation 



-X, ip ^ + Q, t 



T 

^+2 



(5) 



(6) 



leaves the equation ^ invariant and changes the sign 
of V and . Consequently, if Eq.(^) allows for only one 
attractor solution, both quantities will have average value 
zero. Violating |(SJ we loose the symmetry ^ and may 
expect nonzero values for the mean valu.6S_of V and . 
This can be done e.g. by the choice 



E[t) = El cos(ujt) + E2 cos{2ujt + Q) , uj 



27r 

Y 



(7) 



The soliton ratchet effect has been observed in terms 
of the mean soliton velocity V both numerically 0, 0, 
5] and experimentally in an annular Josephson junction 
;8|. However, by differentiating Eq.Q and using Eq.lQJ 
together with the boundary condition it follows that 
[IE El: 



xiptxdx 



(3) 



j/(i) = -^aj\t) ~ QE(t) . 



(8) 
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Thus, for any ac driving force E{t) with zero mean 
the time-averaged value of the total momentum, = 
limt^oo 7 /g J^{T)dT, is zero. Suppose that a kink moves 
with some non-zero mean velocity V. The kink carries 
also some non-zero energy. Hence the progressive kink 
motion should lead to the appearance of an energy flux 
in the system. Thus the system supports a non-zero av- 
eraged energy flow and a vanishing averaged internal en- 
ergy current (total system momentum). We will solve 
the puzzle in this Letter. 

The PDE ^ corresponds to the energy density p 

pMx, t)] = p{x, t) + + cos((^) . (9) 

Using Eq.Q we obtain: 

pt = -ji - a^l -f E{t)ipt ■ (10) 

The last two terms describe energy losses through dis- 
sipation and the energy exchange between the system 
and external driving E{t). However they also contain a 
derivative of an additional current density, as we will 
show below. A central result of this work is that these 
terms describe a new exchange energy current J^{t). 
This exchange current corresponds to an additional en- 
ergy transmission channel provided by a spatially inho- 
mogenous energy exchange between the system, the ex- 
ternal ac driving and the dissipation. Thus, the complete 
current balance equation for the full current J reads: 



J{t)^j'{t) + J^{t) 



(11) 



We consider a large system size, i.e. L ^ Lk, where 
Lk is the kink localization length. Far from the kink 
center the field is oscillating in time while being homo- 
geneous in space. We also assume that the kink trav- 
els only over distances Lp <^ L during one period of 
the ac driving. Because of the external ac force the 
field evolution is uniquely locked to the driver (as it was 
the case for all previous considerations H, Q IE IE 
i.e. assuming that the system allows only for one at- 
tractor). Thus we observe a moving localized excitation 
(kink) which propagates on a background formed by a 
spatially homogenous groundstate. It is convenient to 
separate (/'(x, t) into a localized kink part, (p''{x, t), where 
ip''{x ^ 0,i) = 0; (p''{x L,t) = 2tt, and a back grou nd 
{vacuum) part (p^{t) which depends only on time [lOj: 



ip{x,t)^ip''{x,t)+ip''(t) 



(12) 



The vacuum part alone must satisfy Eq. . Because it is 
also a solution of the system in the absence of a kink when 
Q = 0, it can not contribute to any energy transport j^. 

On the attractor the dynamics of the system is 
given by 

ip\x, t + T)^^\x- VT, t) , (^"(i + T) = ip^'it) , (13) 



where T = 2tt/lij and V is the averaged kink velocity, V = 
(F(t))T = ^ Jo ^i't)dt. Note that all integral system 
characteristics, such as the total energy of the system, 
the kink velocity, and energy currents from Ea. (|ll|l are 
periodic functions of time with the period T. 

Let us compute the full energy current J produced 
by a moving kink. We use the periodicity of the total 
system energy W{t) = p(x, t)dx in time. The amount 
of energy carried through a point x between time t and 
t + T is equal to: 

Aw{x,t)^ I [p{x',t) - p{x',t + T)]dx' ^ 
Jo 

ft+T 

dx' / pt{x',t')dt' . (14) 
Jo Jt 

Here we assume that at x = we have a homogenous 
vacuum state during the full cycle of the ac driving, thus 
there is no energy current through this point. Due to 
the presence of ac driving the system is not time homo- 
geneous and Aw{x,t) in Ea. (|14|l depends on t. With 
@,l|I31) it follows Aw{x,t + T) = Aw{x - VT,t) and 
Aw{x + L,t) = Aw{x, t) and 

J=^^ dt dxj{x,t) , j{x,t) ^ ^Aw{x,t) . {15) 

Using Eas.lf^. ljl^ll and the multiplicative integration rule 
[l3| we obtain from ifTHl . lfTK|l : 



J = V< 



py{x,t)]dx-Lp[^''{t)]h 



(16) 



with (...) — 1/T ...dt. The quantity being averaged 
in the rhs of l|16|) is the difference between the energy 
of system with and without a kink, or simply the kink 
energy W'^ — WQ{t) — W^it), so that the mean total 
current, Ea. p6|l generated by the moving kink reads 



(17) 



It follows that the total energy current has the same sym- 
metry properties as the kink velocity This result 
proves the initial intuitive guess that a moving kink in- 
deed generates a non-zero energy current in the system. 

Because for any ac driving E(t) with zero mean the 
averaged internal current is equal to zero, we arrive 
at the following energy current balance: 



(18) 



In order to obtain an expression for the exchange cur- 
rent density, we use H10|l and (|15|l and arrive at 
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t+T 



f{x,t)^-^J^ dt' dx'^ix',t') , (19) 
(j>{x, t) = = aipt - E{t)ipt , 
J^ = ^^ dtj^ dxj''{x,t). (21) 



(20) 



3 




'0 1 2 3 4 5 6 





FIG. 1: The dependence of the mean exchange current 
on Q for a — 0.2 (sohd hne) and a = 0.05 (dashed hne). 
Circles correspond to the numerical results for J'°*"'. Other 
parameters: £i = £2 = 0.2, uj = 0.1, Q = 27r, L = 500. 

Evidently the exchange current has the same symme- 
try as V, J and . 

We solved the equation numerically in order 
to test the current balance H18|l. A crucial parameter 
is the mesh size h used for the spatial discretization of 
the equation For any finite value of h the internal 
current is nonzero, as also obtained in However 
it scales according to ~ for ft, < 0.1 and vanishes 
in the continuum limit ft — > in full accord with (|18|l . 
We chose ft = 0.1 here, for which the values of J and 
are determined with an error of less than 3%, while 
J^/J ~ 0.003 (for ft = 0.32 used in 5] the latter ratio 
increases to 0.02 for a = 0.2 and 0.18 for a = 0.05). 
The results of numerical calculations of J, Eq. H16(l , and 
, Eq. lj2U, are shown in Fig^ We obtain very good 
agreement with H18|l . Note that L = 500, Lk ~ 10 and 
Lp w 30. According to the collective coordinate the- 
ories 0, Q a- nonzero average kink velocity implies the 
excitation of internal shape modes on the kink. Thus 
the appearance of a nonzero exchange current is another 
manifestation of the excitation of internal kink modes for 
the case of the considered soliton ratchet. 

In Fig|2Ia) we plot the space-time evolution of the 
function (p{x, t) as defined in (|^, which describes the en- 
ergy exchange between the field tp, the external ac drive 
E{t) and the friction term. The energy is exchanged and 
transported in a cyclic way: first the kink absorbs energy 
in its rare tail, then it releases energy in its front, then 
it absorbs energy in its front and finally releases energy 
in the rare tail. In Figl^^b) we plot the space-time evo- 
lution of the energy density p{x,t) ©. The excitation 
of internal shape modes on the kink is clearly observed - 
the kink is much more compressed when moving opposite 
to its average propagation direction as compared to the 
times when it moves in the same direction. 

It is very instructive to apply our approach based on 
the generalized current balance equation Hll|l to the well- 
studied case of a constant force, E = const 0. Due to 
the time homogeneity of the system 1^9(2;, t) = ip{^) , ^ = 
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FIG. 2: Space-time evolution of the soliton ratchet for a = 
0.2 and other parameters as in FigQ (a) Contour plot of the 
function (f){x,t), Eq. I|20|l . (b) Contour plot of the energy 
density J^. 

X — Vt. Still we deal with the motion of a spatially lo- 
calized kink with a finite extension. The internal current 
can be evaluated using Eq.ljSj): 



The internal current density 

j'iO=Vip%^i. (23) 

Similarly the exchange current density can be obtained 
as 

j^(e) = dx{aV\l + VE^,) . (24) 
Jo 

Due to the time homogeneity the kink energy Wk is now 
independent on time, and the total current can be com- 
puted using l(T7|l . We plot the current densities in Fig|31 
The internal current density is single peaked, with the 
peak position corresponding to the position of the kink 
center. The exchange current density ^ 0, which may 
come as a surprise, since in this case the kink moves with- 
out exciting internal shape modes. Yet the kink is a spa- 
tially extended object, and it is this property which leads 
to a nonzero exchange current density. We also com- 
puted the kink velocity, the kink energy, the total current 
and the total internal and exchange currents for the case 
E = 0.2 and a = 0.2, L = 670 and ft = 0.045. The kink 
velocity V = -0.6201, the kink energy Wk = 10.0813 
and thus the total current J = —6.251. The internal cur- 
rent — —6.282, which is off the exact result = — 27r 
(|22|l by an error 5 = 0.001. The total exchange current is 
obtained as = 0.0321. The current balance equation 
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FIG. 3: The densities of the scaled internal current, 0.02j^ (x) 
(thin line), and the exchange current, j^{x) (bold line), for 
the constant force case E — 0.2. The other parameters are 
Q = 27r, L = 500 and a = 0.2. 



appearance of this new exchange current is provided by 
the finite spatial extent of the kink. Even for the case of 
a constant external field, the exchange current is found 
to be nonzero. 

Discretizing space will additionally change the ratio 
between the various currents. For the ac driving case of 
a soliton ratchet, corrections to the internal current inval- 
idate IHIl so that also contributes to the total energy 
flow |5j even for a zero mean drive. While our results are 
important for the general case of spatially extended fields 
coupled to external driving fields or simply other degrees 
of freedom, they also apply directly to the soliton ratchet 
in an annular Josephson junction 

We thank N.R. Quintero for drawing our attention to 
[l^ . and Y. Zolotaryuk and V. Fleurov for helpful dis- 
cussions. 



is satisfied within the same small error S, which is 
an order of magnitude smaller than the computed value 
for . Together with a series of further tests, including 
the analysis of the dependence of the results on L and h, 
we conclude that the exchange current is nonzero for a 
constant field E. It follows from the mismatch between 
the total current J = VWk and the internal current J^. 
Since no internal modes are excited on the kink in this 
case, the appearance of an exchange current (or its den- 
sity) is simply linked to the finite spatial extent of the 
kink, and the corresponding spatially (and thus tempo- 
rally as well) inhomogeneous energy exchange between 
the field if and the external field E and the friction term, 
similar to the soliton ratchet. Notably the exchange cur- 
rent for a constant field E will vanish in the Hamiltonian 
limit E ~ a — Q and also in the overdamped limit a ^ 1 . 

Let us combine both cases from above, i.e. both a con- 
stant and an ac components of the driving E. A careful 
tuning of the parameters leads to an exact cancellation of 
both force components and V — Q. At that point J — 
(|17|l . At the same time, according to Eq. (|22|l . the in- 
ternal current is = —QE'^*-°p / a where £;^*°p is the dc 
component of E. This implies an exact balance between 
the two currents, 

jtotal ^ jE ^ _jl (25) 

It is also possible to have a situation where the sign of the 
total momentum of the system J^, Eq.Q, is opposite to 
the sign of the kink velocity V , and the internal current 
is pumping energy against the kink motion. 

We solved the soliton ratchet puzzle of nonvanishing 
kink velocities and vanishing total momentum (sec also 
[isl and ^3). We identified a new energy pathway which 
is entirely mediated by the spatial and temporal inhomo- 
geneity of the system mediated by the ac driving and the 
spatial extension of the kink. More precisely, internal 
kink modes are crucial in order to sustain this energy ex- 
change current for the soliton ratchet. The reason for the 
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